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Abstract. We study the properties of very faint, sub-L* Lyman break galax- 
ies at z~2-5 — thus far a largely neglected but numerically and energetically 
very important population. We find that the LBG luminosity function under- 
goes luminosity-dependent evolution: the number of luminous galaxies remains 
constant while the number of faint ones grows with time. The total UV luminos- 
ity density increases with cosmic time from at least until reaching a peak 
or a plateau around z~2 — behaviour that is governed by the sub-L* galaxies 
in the LF's "faint tail" . Using broadband SED fitting we find a nearly-linear 
relationship between SFR and galaxy stellar mass at z~2. A typical L* LBG 
at z^2 shows a stellar mass of ~1O 1O M0, remarkably similar to the bimodality 
mass at low redshift. This similarity suggests that the mechanisms responsible 
for the galaxy bimodality at low-z may have also been at play at z^2. 



1. Introduction 

Until very recently, studies of z>l galaxies have focused primarily on luminous, 
vigorously star-forming objects such as submillimetre sources or L~L* LBGs 
that are forming stars at rates of 10s or 100s M Q /yr. Such studies have largely 
neglected the less glamorous but far more numerous faint, sub-L* galaxies that 
are forming stars at much lower rates. Although individually faint, these sub- 
L* objects are very numerous and so collectively they generate more than half 
the total UV luminosity density of the universe at high redshift. They are, 
thus, extremely important contributors to the story of star formation and metal 
enrichment in the early Universe. Moreover, the presence of a break at L* in 
the galaxy LF suggests that galaxies below L* differ substantially from those 
above it. Because of this, our understanding of the mechanisms that drive galaxy 
evolution at high redshift may profit from comparing the properties of the better- 
studied luminous (L > L*) objects with those of the neglected sub-L* galaxies. 
This paper summarizes the results of our recent work aimed at understanding 
these hitherto neglected, individually modest but collectively very important 
objects. They are the long faint tail of the high-z galaxy population. 
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2. Luminosity function and luminosity density 
2.1. Data 

The shape of the galaxy LF bears the imprint of galaxy formation and evolution 
processes. We study the global statistics of sub-L* =3 galaxies at z~2, 3, 4, and 
5 using two deep, large-area, multi-field imaging surveys. 

At z^2, 3, and 4 we use the Keck Deep Fields (KDF) of Sawicki & Thomp- 
son (2005, 2006a, b). The KDF use the very same U n G7ZI filter set and color- 
color selection criteria as are used in the well-known work of lSteidel et al.l ( 19991 . 



2003, 12004 ). However, in contrast to the Steidel et al. work, the KDF reach 



lZii m (AB) =27 , which is 1.5 mag deeper than the Steidel et al. surveys and sig- 
nificantly below L* =3 ; even at z~4 we reach 2 mag fainter than M*. The KDF 
have an area of 169 arcmin 2 divided into three spatially independent regions on 
the sky that allow us to monitor the effects of cosmic variance. The U n GlZI 
filters and color-color selection commonality with the work of Steidel et al. lets 
us confidently select very faint, sub-L* =3 galaxies at z~2, 3, and 4 by relying 
on their extensive spectroscopic characterization of sample selection. It also al- 
lows us to confidently combine our data with theirs, thereby for the first time 
consistently spanning such a large range in galaxy luminosity at high redshift. 

Our z~5 work is bas ed on very d e ep VlrZ 1 im aging of two Subaru Suprime- 
Cam fields presented in llwata et al.l (120061. I2007T). These d ata are wider and 



deeper than earlier work we presented in llwata et al.l ( 20031 ): they now span a 



total of 1,300 arcmin 2 divided over two fields with z^ m (AB) = 26.5 and 25.5. 
The large area covered by these data is essential for determining the abundance 
of the rare, luminous objects, while at the same time the depth of the data 
gives us the necessary ability to probe below L* in the LF. We use VI c z' Lyman 
break selection to select z~5 galaxies and our ongoing spectroscopic program 
has already confirm ed a number of zy5 galaxies, valida ting our z~5 color-color 
selection technique ( Ando et al.ll2004l ; llwata et al.ll2007l ). 



2.2. The differentially evolving luminosity function 

We find that the luminosity function of high-z galaxies undergoes evolution that 
is differential with luminosity. As Fig. 1(a) shows, the bright end of the LF 
remains virtually unchanged, but the faint end evolves from z~5 to at least 
z~3. Our tests show that it is highly unlikely that this observed differential 
evolution is due to some systematic eff ect such as selection bias, our modeling of 
the survey volume, etc. (for details see ISawicki fc Thompson ll2006al ). The effect 



is also quite significant — the statistical probability that the z~4— >3 evolution is 
not differential with luminosity is only 1.5%, and the result is even more robust 
for z~5^3. We conclude that the differential, luminosity- dependent evolution is 
very likely real. If so, it must reflect some real, intrinsic evolutionary differences 
between faint and luminous LBGs. Understanding what causes these differences 
will give us insights into what drives galaxy evolution at high redshift. 

2.3. Luminosity and star formation rate densities 

The importance of the LF's faint end extends into measurements of the cosmic 
luminosity density and the cosmic star formation rate that is often derived from 
it. For reasonable faint-end slopes of the LF the bulk of the luminosity density 
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Figure 1 . (a) Luminosity-de pendent evolution o f the UV luminosity func- 
tion. The z~5 results are from llwata et al.1 (l2007[) while the z ~4, 3, and 2 
points are f rom the KDF work o f lSawicki k. Thompson I (|2006al ) and include 
the data of ISteidel et~aTl (|1999l ) at the bright end. The z~2.2 LF may be 
subject to a systematic bias that could be suppressing the number density 
shown by up to a factor of tw o; however, the z~3 and z~4 points are ex- 
tremely robust, as discussed in ISawicki &: Thompson I (|2006af h (b) The UV 
luminosity density of the Universe and the corresponding star formation rate 
density (uncorrected for dust). Note that the divisions by luminosity are de- 
fined on the basis of the non-evolving z^3 L^ =3 . The z~5, 4, 3, and 2 points 
are from our work, while the lower- z data are from GALEX. The bulk of the 
UV luminosity density at high redshift, z>2, comes from the faint, thus far 
largely neglected galaxies in the range L=(0.1-1)L* =3 . 



resides in galaxies that are fainter than L* and it is only with deep, multi- 
field, large-area surveys such as ours that it is possible to reliably measure the 
contribution of this dominant population. 

Figure 1( b) shows the evolution of the UV lum inosity density using data 
from GALEX (|Arnouts et al.ll2005l: IWyder et al.ll2005l ). Keck Sawicki & Thomp- 
son (2006b), and Subaru ( Iwata et al. 20071 ). The total UV luminosity density, 
i.e., luminosity density due to galaxies of all luminosities, rises from early epochs, 
z>5, experiences either a plateau or a broad peak at z~3-l, and then a decline 
to z=0. This behavior of the total luminosity density is dominated not by lumi- 
nous galaxies, but by sub-L* =3 objects. Indeed, it is galaxies within the rather 
narrow luminosity range L=(0.1-1)L* =3 that dominate at high redshift, z~2-5, 
and it is the rapid evolution of this population that drives the evolution in the 
total UV luminosity density. Thus far, this faint population has been largely 
neglected in high-z follow-up studies which have focused on L > L* z= ^ LBGs. 
In contrast to z>2, at lower redshifts, z<l, the total UV luminosity density is 
dominated by very faint galaxies with L<0.1L* =3 . The L=(0.1-1)L* =3 galax- 
ies that dominated at z>2 are still important but no longer dominant. This 
switch in the luminosity of the galaxies that dominate the total luminosity den- 
sity is reminiscent of galaxy downsizing, although the effect we see so far is a 
downsizing in UV luminosity rather than galaxy mass. 
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3. Physical properties from SED fitting of sub-L* galaxies at z~2 



3.1. The galaxy sample and SED fitting technique 

As we have discussed above, sub-L* LBGs dominate over their more luminous 
cousins, but evolve differently from them as a population. We now wish to 
know what makes them physically different. We attempt to address this issue 
through broadband spectral energy distribution (SED) studies. Unlike in the 
LF work, here we do not need large samples of objects but instead must have 
deep multiwavelength photometry that spans from the UV to beyond the Balmer 
break in the rest frame. We need to go particularly deep because we wish to 
study galaxies significantly fainter than L* at z>2. No spectroscopic samples of 
such very faint LBGs currently exist and so we must use photometric redshifts. 

Our first foray into this area is at z~2 and uses the HST WFPC2+NICMOS 
^300^450^606^814^110^160 data of the northern Hubble Deep Field (HDF). These 
data are ideal for us as they are very deep, include the t7-band imaging essential 
for selecting z~2 LBGs, and span the age-sensitive Balmer break at the target 
epoch. 

After smoo thing the multiband im ages to a common PSF and obtain- 
ing SExtr actor ( Bertin &; Arnouti 19961 ) photometry, we constructed a [/-band 
drop-out sample using the LBG color-color selection criteria defined by Stei- 
del et al, (1996) for the HDF. We then applied a photometric redshift cut, 
1.8<Zp/j t<2.6, where the photometric redshifts are determined as part of the 
SED-fitting procedure (see below). The lower limit of this photo- z cut eliminates 
potential low-z interlopers and the higher limit excludes high-z LBGs for which 
the C/3oo-B45oV606-^8i4Lno-£fi6o fi^er se ^ does not span the crucial, age-sensitive 
Balmer break. This procedure (two-color selection followed by a photo-z cut) 
is very similar to the procedure (two-color selection followed by spectroscopy) 
used bv lSteidel et al.l in their work. Our resulting sample of ~70 objects has a 
m ean redshift Z p/ w=2 .3 an d in other ways also c losely mimics the BX samples 
of ISteidel et all ( 20031 ) and IShaplev et all (|200 5). However, it contains much 
fainter galaxies, reaching down to R(AB)=27^ 

Our SED fitting p rocedure follows the no w well-established approach first 
developed f or LBGs by Sawicki &: Yee I (I1998T) a nd su bsequ e ntly used by many 
others (e.g.. iPapovich et al.ll200ll : IShaplev et al. | [200ll . 120051 : llwata et al.ll2005T ). 
We compa re the observed -B^ o Vftoftls uJiioffifirEl galaxy photometry with pre- 
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dictions of Bruzual Charlotl (|2003l ) spectral synthesis models attenuated by 



Hboo) dust. We fix metallicity and star formation history and 
fit for burst age, star formation rate, stellar mass, dust reddening, and redshift. 
While for brighter samples redshifts are usually constrained from spectroscopy, 
here we must keep redshift as a free par ameter. Howe v er, our tests as well as 
previous experience with this approach ( Sawickil 120021 : lHall et aL 200 ll ) show 
that lack of spectroscopic redshifts does not significantly affect the results. 



1 As in ISteidel etHl <|1996fl . the 7? magnitude is defined via averaging the V(jo6 and /gi4 fluxes. 

2 While C/300 is essential for selecting our LBG sample, we omit it from the SED fitting to 
avoid contaminating our SED analysis with the stochastic effects of intergalactic hydrogen 
absorption. 
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Figure 2. (a) Stellar mass as a function of UV m agnitude at z~2.3. Results 
for brighter galaxies fitted bv lShaplev et alj (|2005l ) are sh own along with our 
sampl e. Both samples were fitte d with constant SFR Bruz ual fc Chariot 
(|2003l ) models attenuated with iCalzetti et all (|2000| ) dust, (b) Stellar mass 
as a function of (dust-corrected) star formation rate. In both panels solid 
lines show fits to the HDF data. 



3.2. A mass-SFR relation and a characteristic mass 

Figure 2a shows that the bulk of the galaxies in our HDF sample follows a re- 
lation between stellar mass and UV l umin osity and this relation extends into 
the brighter galaxies of lShaplev et al.1 ( 20051 ). There appears to be good correla- 



tions between stellar mass and UV luminosity (Fig. 2a) and stellar mass and SFR 
(Fig. 2b). Outliers in both samples turn out to be galaxies with constant-SFR 
best-fit ages apparently older than the age of the Universe at z=2.3 (trian- 
gles) and galaxies with very young ages age/it <30Gyr (squares). However, most 
objects show a fairly tight correlation which can be described by the relation 
log(M stars /M ) = 9.0 + 0.861og[S"F J R/(M Q yr- 1 )] (solid line in Fig. 2b). The 
inclusion of the outliers in the fit does not significantly affect this result. 

The L* knee in the z~2.2 LF (M 1 * 70M =-20.6 at z~2.2: Sawicki & Thomp- 
son, 2006a) corresponds to a stellar mass M<j tar . s (L*)~lO lo M (Fig 2a). This 
z~2 characteristic mass is remarkably similar to the 3x10 10 Mq characteristic 
mass that at low reds hift marks the transition between the two distinct galaxy 
'bimodality" families ( Kauffmann et aTl 12003]) . This similarity of cha racteristic 



masses suggests that the processes (see, e.g.. |Pekel fc Birnboim|[20 06) responsi 



ble for the observed galaxy bimodality at low z may have also been operating 
at z~2. Moreover, the near-unity slope of the logarithmic relation between stel- 
lar mass and UV luminosity (Fig. 2a) suggests that the bulk of the detectable 
stellar mass in z~2.2 LBGs resides in the thus-far poorly studied galaxies with 
L*>L>0.1L* , just as does the bulk of the total UV luminosity density. 

The outliers with extremely young fit ages (square symbols) are apparently 
either undermassive for their SFRs or overluminous for their masses. In the 
former case they may move up onto the mass-SFR relation as they age; in the 
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latter they could be experiencing strong star-forming bursts that dramatically 
elevate their normally low luminosities. Outliers with very old best-fit ages 
(triangles), are also consistent with variable star formati on histories: constant 
star formation histor ies tend to yield the largest ages (e.g. , ISawicki & Yeell998l : 
iPapovich et al. 200ll ) and thus the old ages here may simply be telling us that 
star formation in these galaxies was more intense in the past. All these scenarios 
advocate that at least some z~2 LBGs, both faint and luminous, have variable 
star formation histori e s — c onsistent with one of the scenarios proposed by 
Sawicki Thompson (2006a) to explain the evolving LBG luminosity function. 



4. Concluding remarks 

The thus-far largely unexplored sub-L* LBGs dominate the UV luminosity den- 
sity of the universe at z>2. Moreover, their stellar masses appear to correlate 
with UV luminosities down to very faint magnitudes at least at z~2 and so, 
just as they dominate the UV luminosity density, these objects may also contain 
a substantial fraction of the formed stellar mass by that redshift. At the same 
time, however, the differentially-evolving luminosity function suggests that these 
sub-L* LBGs differ substantially from their more luminous cousins. It is not yet 
clear what mechanism is behind this differential evolution, but the bursty na- 
ture of at least some of the z~2 galaxies means that changes in the timescales 
of star-burst episodes may play a role. We will address this and other issues as 
we extend our SED-fitting and other differential studies to higher redshifts. 
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